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Abstract. In this paper we propose a copula contagion mixture model for correlated default 
times. The model includes the well known factor, copula, and contagion models as its special 
cases. The key advantage of such a model is that we can study the interaction of different 
models and their pricing impact. Specifically, we model the marginal default times to follow 
some contagion intensity processes coupled with copula dependence structure. We apply 
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the total hazard construction method to generate ordered default times and numerically 
compare the pricing impact of different models on basket CDSs and CDOs in the presence 
of exponential decay and counterparty risk. 
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^ ! 1 Introduction 

O _ 

The recent financial crisis has profound impact on the financial systems in the US, UK, and 
other major markets. Some giant banks and insurance companies either collapsed or had to 
be bailed out by the national governments. The excessive risk exposure of many banks to 
collateral debt obligations (CDOs) and credit default swaps (CDSs) has played the key role 
in this financial crisis. One may list many causes which have attributed to and aggravated 
the crisis, however, in this paper we only focus on the impact of correlation modelling on the 
pricing of these portfolio credit derivatives. 

CDOs and CDSs had phenomenal growth in recent years until this financial crisis. The 
key in pricing and hedging these portfolio credit derivatives is the correlation modelling. 
There are mainly three approaches in the literature: conditional independence, copula, and 
contagion. Factor models are most popular due to their semi-analytic tractability. Many 
effective algorithms have been developed to characterize the portfolio loss distribution, see 
Andersen et al. (2003), Hull and White (2004) for recursive exact methods, and Glasserman 
(2004), Zheng (2006) for analytic approximation methods. Factor models may underestimate 
the portfolio tail risk and economic capital, see Das et al. (2007). Copula models are also 
popular, especially the Gaussian copula which is used in CreditMetrics, see Li (1999). Some 
copulas (Archemedian and exponential) are good to model extreme tail events and simulta- 
neous defaults. There are some active recent debates on the usefulness of copulas in financial 
modelling and risk management, see Mikosch (2006) and many discussion papers on the same 
issue. Contagion models study the direct interaction of names in which the default intensity 
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of one name may change upon defaults of other names and "infectious defaults" may de- 
velop, see Davis and Lo (2001), Jarrow and Yu (2001), Yu (2007). It is in general difficult to 
characterize the joint distribution of default times due to the looping dependence structure. 
For homogeneous portfolios there is a closed form formula for the density function of ordered 
default times, see Zheng and Jiang (2009). Monte Carlo method is often used to price CDOs 
and basket CDSs no matter which correlation model is used and provides benchmark results 
to test efficiency and accuracy of analytic and numerical algorithms. 

It is interesting to know which model one should choose in pricing portfolio credit deriva- 
tives. We know different models give different values. If one uses the Gaussian copula, the 
swap rate for senior tranche of CDO is low due to the thin tail distribution of portfolio loss, 
on the other hand, if one uses the contagion model, the swap rate for the same senior tranche 
is much higher. However, one cannot simply say the contagion model is preferable to the 
Gaussian copula because it provides higher swap rates for senior tranche. It all depends 
on the underlying model assumptions. These correlation models are defined under different 
frameworks and are difficult to compare directly their pricing impact. It is therefore benefi- 
cial to have a unified model which covers all three known models as special cases. One may 
then extract the information of the interaction of these models and may give a more balanced 
view on which model one should choose for a specific application. 

In this paper we suggest a general copula contagion mixture model which includes factor, 
copula, and contagion models as its special cases. The key advantage of such a model is 
that we can study the interaction of different models and their pricing impact. Specifically, 
we model the marginal default times to follow some exponential decay contagion intensity 
processes coupled with some copula dependence structure. This is not a Markov process 
model and cannot be solved with the standard Kolmogorov equations or matrix exponen- 
tials, see Frey and Backhaus (2004), Herbertsson and Rootzen (2006). Although there are 
analytic pricing formulas for some special cases, we choose to use the Monte Carlo method to 
price CDOs and basket CDSs, which is reliable, accurate and efficient with some optimized 
numerical procedure. 

The paper is organized as follows: section 2 describes the copula contagion mixture model 
and the relation with the known models, section 3 applies the model to price CDOs and basket 
CDSs and discuss the impact of interaction of different models, section 4 concludes. 

2 The Model 

Let (f2, T, {Tt}t>o, P) be a filtered probability space, where P is the martingale measure and 
{•7-"}i>o is the filtration satisfying the usual conditions. Let Tj be the default time of name i, 
Ni(t) = l{ Ti <t} the default indicator process of name i, T\ = a(Ni(s) : s < t) the filtration 
generated by default process TVj, i = 1, . . . , n, and T% = T\ V ... V the smallest a-algebra 
needed to support n, . . . , r n . Assume that Tj possesses a nonnegative Ft predictable intensity 
process \{t) satisfying E[Jq X(s)ds] < oo for all t, and the compensated process 

rtATi 

Mi(t) = Ni(t) - / Xi(s)ds 
J o 

is an Ft martingale. Given Tj = tj, j G Jk = {ji, ■ ■ ■ ,jk} C {1, . . . , n}, satisfying = tj < 
tj 1 < . . . < tj k and Tj > t > tj k for i J^, the conditional hazard rate of Tj at time t is given 
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by | 

W\tj k ) = |m -^P(t <n<t + At\ Tj = tj,j e J fc ) 

where t j h is a short form for (tj 1 , . . . , tj k ) . 

In copula modelling of default times it is normally assumed that intensity processes are 
independent of default times of other names, i.e., \i(t\tj k ) = Xi(t) for all t. The marginal dis- 
tribution functions of default times r% are given by = P(t{ < t) = E[l— exp(— J Q * Aj(s)ds)] 
(if Xi are stochastic processes) and standard uniform variables i^(rj), z = l,...,n, have a 
joint distribution function C, a given copula. It is easy to generate default times Tj with 
the Monte Carlo method. One can simply first generate standard uniform variables U with 
copula C and sample paths of Aj and then find the default times Tj by 



inf jt > : y Aj(s)ds > j 



where £7j = ln(l — Ui), i = 1, . . . , n, are correlated standard exponential variables. In partic- 
ular, if Xi(t) = di, then Tj = -Ej/dj. 

There has been extensive research in literature on factor modelling of default times. These 
models are all special cases of copula modelling of default times. For example, the well-known 
Gaussian factor model is given by 

Xi = pZ + - P 2 Zi 

where Z, Z\, . . . ,Z n are independent standard normal variables and p is a constant satisfying 
\p\ < 1. Z is often interpreted as a systematic factor and Zi idiosyncratic factors. If we set 
U = $(Xi), where $ is the standard normal distribution function, then the distribution of 
U is a Gaussian copula given by 

C(m, ...,u n ) = $ mi r(<I>"Vi)> . . . , $ _1 K)) 

where $ m ,r is the n-variate standard normal distribution function with mean vector m = 
and correlation matrix F that has diagonal elements 1 and all other elements p 2 . Factor 
models are appealing from model interpretation and conditional independence point of view. 
The corresponding copulas may have some complex forms, but in general it is easy to generate 
correlated standard uniform variables due to the special structure of factor models. From 
mathematics point of view there is no need to treat them separately if we know how to 
generate standard uniform variables U from given copulas C . 

In contagion modelling of default times the intensity processes Xi(t\tj k ) depend on default 
times of other names and standard uniform variables Fj(rj) are assumed to be independent. 
The marginal distribution functions Fj of default times Tj cannot be simply expressed in terms 
of Aj(s) as information of other default times is needed to characterize the whole intensity 
process paths and there is a "looping" phenomenon. Although it is difficult to characterize 
the marginal and joint distributions of default times it is easy and straightforward to generate 
ordered default times t 1 with the total hazard construction method. One can first generate 
independent standard uniform variables Ui and set Ei = — ln(l — U), i = 1, . . . , n, then find 
default times one by one as follows: To find the first default time r 1 and the corresponding 
name j\ , set 

h = argmin j=lj n |tj- > : X j (s)ds > E^ 



and r 1 = tj 1 and J\ = {ji}, where Xj(s) are unconditional hazard rates of names j at time 
s. To find the kth default time r k and the corresponding name j& for k > 2, set 

j fc = argmin^ Jfe i |^ > r ifc _ 1 : jf X^tj^da > Ej\ 

and r fe = t jk and J fe = J k -i U {j fc }. 

We suggest a copula contagion mixture model which covers both copula model and con- 
tagion model as special cases. Specifically, we assume that the intensity processes Aj may 
depend on default times of other names and standard uniform variables i^(rj) have a joint 
distribution C. This is a natural generalization of pure copula models and pure contagion 
models. One can easily generate default times with the total hazard construction method. 
The only difference with the pure contagion model is that one generates standard uniform 
variables Ui from a given copula C, not necessarily from the product copula which corre- 
sponds to the pure contagion model. The key advantage of this new mixture model is that, 
instead of studying three well known models in isolation, we can explore their interaction and 
their joint pricing impact on CDOs and basket CDSs. 

We now impose some structure to the intensity processes. To simplify the notation and 
highlight the key point, we assume a homogeneous portfolio. The discussion is the same 
for general heterogeneous intensity processes except the expression is more complicated. We 
assume the intensity processes have the following structure 

A i (t) = a(l+ fl ce ~ d{t ~ T]h {T 3 <t} I , i = l,-..,n, (1) 

\ i=^3+i J 

where a, c, d are positive constants. (These parameters can be deterministic functions of 
t or even some stochastic processes, the discussion is essentially the same, see Zheng and 
Jiang (2009).) a is the unconditional default intensity, c is the contagion rate, and d is the 
exponential decay rate. When d = we may introduce the default state space and use the 
Markov Chain to study the joint distribution of default times. Apart from this extreme case 
the intensity processes ([1]) are non-Markov. 

For homogeneous intensity processes ([1]) without exponential decay (d = 0) we can sim- 
plify the total hazard construction method. This is because we only need to know the number 
of defaults at time t but not the identities of names which have defaulted. We can generate 
T k as follows: 

Step 1. Generate correlated standard uniform variables Ui, i = 1, . . . , n, from the copula C. 
Step 2. Set Ei = — ln(l — Ui), i = 1, . . . , n, and sort Ei in increasing order to get E* with 
E%<E%<...< E*. 

Step 3. Find ordered default times r fc by setting 

TP* C 1 * P* 

= r k =r^+ ,f* fOT k = 2,...,n. (2) 

a a(l + (k — l)c) 

The density function of the feth default time T k is given by (for d = 0) 

fc-i 

f Tk {t) = Y. a ^ ae ~ p ° at 

3=0 
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where (3j = (n — j)(l + jc) and ctkj are constants depending on c,k,j,n and have explicit 
expressions. For example, f T i (t) = nae~ nat which shows that the contagion has no influence 
on the first default time, and 



n (t+S-n)t (- e ~ nat + e-("- 1 )( 1 + c )^), if c ^ l/(n - 1) 
(na) 2 te- nat , ifc = l/(ra-l) 



which implies that the contagion affects the second and all subsequent default times. We can 
then derive the analytic pricing formulas for basket CDSs and CDOs, see Zheng and Jiang 
(2009). 

For general homogeneous intensity processes (pQ) we cannot use ([2]) to generate ordered 
default times. The computation is slightly more involved. Steps 1 to 2 are the same and 
so is the first default time r 1 . Assume ordered default times r 1 , . . . ,r k ~ l have already been 
generated for some k > 2. Now we want to generate T k . Let r k ~ 1 <t<T k . The total hazard 
accumulated by name k at time t is 



fc— 1 ,.-ri / 7 — 1 \ ,.t / fe— 1 



X k {s)ds = J2j a^l + Y,ce- d( - s - T,) jds + J a ^1 + ^ ce"^"^ J ds 
where r° = and Y^=i = by convention. Simplifying the above expression we get 

J\ k (s)ds = at + - e~ d ^) . 

k 

The T k is determined by the relation J Q T Xk(s)ds = E k . Define 



k-l 

ac 

i=i 

Then r k is a root of nonlinear equation Ff.(t) = 0. Since F!(t) > and F','(t) < function F k 
is strictly increasing and strictly concave. Observe also that from i ? / c _i(r fc_1 ) = we have 

F fc (T^) = ar k ^ + f E (! - e^"^) ~ K = E^ x - E* h < 

i=i 

and Fk(oc) = oo. There is a unique root of equation F k (t) = on the interval [r fe_1 ,oo). 
The special structure of function F k guarantees that the Newton algorithm with an initial 
iterating point r fc_1 converges quadratically to the root r k . We can now summarize Sept 3 
in the presence of exponential decay rate d > as follows. 

Sept 3 ' . Set r 1 = E\ja and find the kth default time T k by solving numerically the equation 
Fk(t) = with the Newton algorithm and the initial iterating point r k ~ 1 for k = 2, . . . , n. 

We now discuss the impact of exponential decay rate d on ordered default times T k . From 
F' k (t) = a + ac 

g-rf(t-r') w know that 

is a strictly decreasing function of d for 
t > T k ~ 1 . If d = we have F!(t) = a + ac(k — 1) and F k is a linear function 

F k (t) = E* k _ x -E* k + (a + ac(k - l))(i - r*" 1 ). 



The kth default time r k is given by ([2]) as expected. If d = oo we have F' k (t) = a and F k is 
again a linear function 

F k {t) = El^-El+a(t-r k - 1 ). 

The fcth default time is given by r k = r fc_1 + (E k — E^_^)/a, or equivalently, r k = E%/a, 
which corresponds to the case when there is no contagion effect. For any other d the kth 
default time r k lies between these two extreme cases. We conclude that the smaller the 
exponential decay rate, the stronger the contagion effect and the sooner the ordered default 
times, which makes CDO and basket CDS riskier and demands higher spreads. 



3 Numerical Tests 

We can now value the basket CDS and CDO with the copula contagion mixture model. 
We assume homogeneous intensity processes ([1]) to simplify the computation, but the same 
method can be applied to general intensity processes. For both basket CDS and CDO we 
assume that T is the maturity of the contract, t\ < t<i . . . < tjy are swap rate payment dates, 
to = is the initial time and tjy = T is the terminal time, R is the recovery rate, r is the 
riskless interest rate, and B{t) = e~ rt is the discount factor at time t. 

To price basket CDS we assume S k is the annualized kth. default swap rate. The expected 
value of the contingent leg at time is equal to 

E[(l-R)B(r k )l {Tk < T} ] 

and that of the fee leg with accrued interest is equal to 

N 



S k E 



i=l 



£ (fe - ti_l)B&)l {Tfc>ti} + (T k ~ t t ^)B(T k )l {u _ 1<n 



k <U} 



We can easily find the swap rate S k with the Monte Carlo method by generating ordered 
default times r k . 

To price CDO we assume ki, I = 0, . . . , M — 1, are attachment points of tranches I with 
= ko < ki < . . . < ku = 1 5 = k\ — k\^\ are tranche sizes for I = 1,...,M, the 
cumulative percentage portfolio loss at time t is given by 



- 1 } 



fc=i 



with r° = and r n+1 = oo, the cumulative tranche I loss at time t is given by 

L i{ 1 ) = (Hi) ~ ^-l) 1 {fei_ 1 <i(*)<fc ! } + AM{L(t)>fc ; }- 

Assume Si is the swap rate of tranche I. The expected value of the contingent leg for tranche 
I loss at time is given by (note L/(0) = 0) 



E 



N 



^B(U){Li(ti)-Li(U-i)) 



p/c 


tranche 


0.0 0.3 3.0 


0.0 


0-0.15 
0.15-0.3 
0.3-1 


0.0740 0.0890 0.2360 
0.0000 0.0003 0.1052 
0.0000 0.0000 0.0199 


0.5 


0-0.15 
0.15-0.3 
0.3-1 


0.0682 0.0843 0.1553 
0.0042 0.0164 0.1020 
0.0001 0.0022 0.0596 


0.9 


0-0.15 
0.15-0.3 
0.3-1 


0.0326 0.0373 0.0488 
0.0147 0.0242 0.0439 
0.0044 0.0157 0.0405 



Table 1: CDO rates with the Gaussian copula contagion mixture model. 



and that of the fee leg for tranche I is 

N 



J2(U - U-i)B(U)(Aki - L,(t0) 



i=i 



We can again easily find the swap rate Si with the Monte Carlo method. 

To generate ordered default times we must first generate correlated standard uniform 
variables Ui, i = 1, . . . , n. We use three different copulas to generate f/j. The first one is the 
product copula and Ui are simply independent standard uniform variables. The second one 
is the exponential copula and XJ{ are generated as follows: first generate n + 1 independent 
exponential variables To, T± . . . , T n , where To has parameter Co and T±, . . . , T n have parameter 
ci, then set Si = min(To, Tj), and finally define £/j = 1 — exp(— (co + cfiSi), i = 1, . . . , n. This 
is the simplest exponential copula which models simultaneous jumps as well as individual 
jumps, see Giesecke (2003) for more details on exponential copulas and Xu and Zheng (2009) 
for their applications in modelling portfolio asset price processes. The third model is the 
Gaussian copula and U are generated as follows: first generate n + 1 independent standard 
normal variables Z, Z±, . . . , Z n , then set Xj = pZ+ y/l — p 2 Zi, i = 1, . . . , n, and finally define 
Ui = &(Xi), i = 1, . . . ,n. This is the most popular model used in financial institutions for 
pricing portfolio derivatives. 

We have used the following data in numerical tests: number of names n = 40, riskless 
interest rate r = 0.05, time to maturity T = 3, number of payments N = 6 with equally spaced 
time intervals, unconditional intensity rate a = 0.01, recovery rate R = 0.5, exponential decay 
rate d = 0, and number of simulations is 1 million. 

Table 1 lists CDO rates computed with the Gaussian copula contagion mixture model 
with different c and p. We can see that swap rates increase if c increases, which is expected 
as higher c causes higher contagion and more defaults, c = corresponds to the Gaussian 
factor model. As p increases swap rates for equity tranche decrease while those for mezzanine 
and senior tranches increase, a well known fact, p = corresponds to the pure contagion 
model (or the product copula contagion mixture model) and we see c has huge impact on 
swap rates for mezzanine and senior tranches. When both c and p are positive, we see swap 
rates for senior tranche are greater than those with the pure contagion model (p = 0) and 
the pure factor model (c = 0). It is interesting to note that swap rates for mezzanine tranche 
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0.3 


3 


k 


ProdC 


ExpC 


GausC 


ProdC 


ExpC 


GausC 


ProdC 


ExpC 


GausC 


1 


0.2024 


0.1575 


0.1153 


0.2024 


0.1575 


0.1153 


0.2024 


0.1575 


0.1153 


2 


0.0634 


0.0697 


0.0508 


0.0769 


0.0811 


0.0573 


0.1401 


0.1249 


0.0855 


5 


0.0010 


0.0026 


0.0105 


0.0052 


0.0104 


0.0197 


0.0836 


0.0866 


0.0620 


10 


0.0000 


0.0000 


0.0014 


0.0000 


0.0001 


0.0072 


0.0486 


0.0582 


0.0492 


20 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0016 


0.0163 


0.0263 


0.0369 


30 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0003 


0.0024 


0.0061 


0.0274 


0-0.15 


0.0740 


0.0742 


0.0682 


0.0890 


0.0923 


0.0843 


0.2360 


0.2218 


0.1553 


0.15-0.3 


0.0000 


0.0000 


0.0042 


0.0003 


0.0011 


0.0164 


0.1052 


0.1246 


0.1020 


0.3-1 


0.0000 


0.0000 


0.0001 


0.0000 


0.0000 


0.0022 


0.0199 


0.0314 


0.0596 



Table 2: Comparison of basket CDS and CDO rates with the product copula contagion 
mixture model (ProdC), the exponential copula contagion mixture model (ExpC) with cq = 
0.01 and c\ = 0.1, and the Gaussian copula contagion mixture model (GausC) with p = 0.5. 

decrease as p increases when c = 3, an opposite phenomenon to the case when c = 0. This 
is not surprising because when c = 3 the default intensity increases quickly for surviving 
names and many more names are likely to default, in other words, the mezzanine tranche 
behaves increasingly like the equity tranche, and therefore as p increases the corresponding 
swap rates actually decrease. For the same reason the senior tranche behaves increasingly 
like the mezzanine tranche and its swap rates increase and then decrease as p increases. We 
have also done numerical tests for a = 0.1 and found that all tranches behave like the equity 
tranche and swap rates decrease as p increases even when c = 0.3 

Table 1 may shed some light on the cause of recent financial crisis. Before the full scale 
credit crunch, the housing and stock markets were booming, the credit was cheaply and easily 
available, few individuals and companies defaulted, and default rates from rating agencies 
were very low. Portfolio credit derivatives such as synthetic CDOs were in high demand. The 
Gaussian factor model (corresponding to c = in the table) was the most popular model used 
in financial institutions to price these securities. Table 1 shows that for the senior tranche 
(0.3-1) the risk is almost negligible for p = 0.5, and is still very small even for unlikely p = 0.9. 
It seemed that the underwriter of CDS for CDO senior tranches could make huge profit from 
premium fees with little risk, almost like "free lunch with vanishing risk". However, when 
there is contagion, which is the case for synthetic CDOs (the actual loss can be many times 
over the nominal loss), the risk for the senior tranche is much higher even in good economy 
(a = 0.01). When c = 3 and p = 0.5, the swap rate for the senior tranche is 0.0596, in 
sharp contrast to 0.0001 when c = and p = 0.5. Mis-pricing of synthetic CDOs could be 
one of the causes which led to the financial crisis of underwriters of CDS for these synthetic 
portfolio credit derivatives. 

Table 2 lists basket CDS and CDO rates computed with three copula contagion mixture 
models. The copulas used are the product copula, the exponential copula with cq = 0.01 and 
ci =0.1 (individual jumps is much more likely than a systematic jump), and the Gaussian 
copula with p = 0.5. The results for basket CDS are mixed with no single model dominating 
the others in pricing. Contagion has no influence to the 1st default CDS rate and the product 
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A; 



d = d = l d = 10 d = 100 d = oo 



1 

2 
5 

10 
20 
30 



0.1153 0.1153 0.1153 0.1153 0.1153 

0.0855 0.0761 0.0564 0.0514 0.0508 

0.0620 0.0482 0.0175 0.0111 0.0105 

0.0492 0.0348 0.0053 0.0017 0.0014 

0.0369 0.0230 0.0008 0.0001 0.0000 

0.0274 0.0137 0.0001 0.0000 0.0000 



0-0.15 
0.15-0.3 
0.3-1 



0.1553 0.1323 0.0810 0.0696 0.0682 
0.1020 0.0727 0.0127 0.0048 0.0042 
0.0596 0.0328 0.0012 0.0002 0.0001 



Table 3: Comparison of basket CDS and CDO rates with exponential decay Gaussian copula 
contagion mixture model, data used are a = 0.01, p = 0.5, and c = 3. 

copula produces the highest rate. When there is no contagion (c = 0) or low contagion 
(c = 0.3) the Gaussian copula dominates the swap rates for all k but the first few. When there 
is high contagion (c = 3) the results are more mixed with the Gaussian copula dominating 
for large k and the other copulas for small k. The results for CDO are also mixed. For equity 
tranche the product and exponential copulas produce similar rates which are higher than 
those from the Gaussian copula. For mezzanine and senior tranches the Gaussian copula 
gives much higher rates than the other two copulas do except when contagion is high (c = 3) 
and the rates from the other two copulas are also increased significantly. The difference 
between swap rates using different copula contagion mixture models is substantial. 

Table 3 lists basket CDS and CDO rates with the exponential decay Gaussian copula 
contagion mixture model. The data used are a = 0.01, c = 3, p = 0.5, and different decay 
rates, d = corresponds to the Gaussian copula contagion mixture model without decay and 
d = co to the case without contagion effect. It is clear that as d increases, basket CDS and 
CDO rates decrease. The exponential decay has much greater impact to the fcth default rates 
for larger k than for smaller k. The same phenomenon is observed for CDO rates, that is, 
the exponential decay has much greater impact to senior tranche rates than to junior tranche 
rates. Basket CDS and CDO rates are highly sensitive to exponential decay rates d, which 
requires an accurate estimation of d in calibration if one is to use it in pricing. 

We have done some numerical tests on pricing of CDOs and basket CDSs in the presence 
of counterparty risk. Assume the intensity processes (pQ) for underlying names in the portfolio 
and assume the intensity process of default time t b of the counterparty is given by 



where as is the unconditional default intensity and cb the contagion rate. Note that the 
hazard rate process Xb of the counterparty is influenced by defaults of names in the portfolio, 
but not vice versa. This follows the observation in Leung and Kwok (2005) and Yu (2007) 
that the contagion of the counterparty on underlying names does not affect CDS pricing. 
To price basket CDS we only need to compute the expected value of the contingent leg 
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c 


0.0 


0.3 


3.0 


p 


tranche 


GausC 


GausCCR 


GausC 


GausCCR 


GausC 


GausCCR 


0.0 


0-0.15 


0.0740 


0.0740 


0.0890 


0.0889 


0.2360 


0.2347 




0.15-0.3 


0.0000 


0.0000 


0.0003 


0.0003 


0.1052 


0.1027 




0.3-1 


0.0000 


0.0000 


0.0000 


0.0000 


0.0199 


0.0188 


0.5 


0-0.15 


0.0682 


0.0680 


0.0843 


0.0841 


0.1553 


0.1521 




0.15-0.3 


0.0042 


0.0040 


0.0164 


0.0160 


0.1020 


0.0968 




0.3-1 


0.0001 


0.0001 


0.0022 


0.0020 


0.0596 


0.0500 


0.9 


0-0.15 


0.0326 


0.0326 


0.0373 


0.0364 


0.0488 


0.0421 




0.15-0.3 


0.0147 


0.0144 


0.0242 


0.0232 


0.0439 


0.0355 




0.3-1 


0.0044 


0.0040 


0.0157 


0.0137 


0.0405 


0.0291 



Table 4: CDO rates with the Gaussian copula contagion mixture model (GausC) and the 
same model with counterparty risk (GausCCR). 



and the fee leg at time 0, given respectively by 



E 



{l-R)B(T k )l {Tk ^ rB > Tk} 



and 



' N 

£((*, 



)B(f,)l( 



T k >ti,T- 



>>ti} + ~ i i-l) B ( T<: ) 1 ft-i<rk<i„rfl>r*}) 



Similarly to price CDO tranche I we only need to compute the expected value of the 
contingent leg and the fee leg at time 0, given respectively by 



E 



N 



^£(i;)(L^)-M**-l))V> ti } 



i=l 



and 



SiE 



N 



J2(U - U-i)B(U)(Aki - L,(ti))l {T B >ti} 



We can easily find swap rates with the total hazard construction method by generating 
ordered default times T k and counterparty default time t b . 

Table 4 lists the swap rates of all tranches with and without counterparty risk. We have 
used the data as = a/10 = 0.001 and cb = c. The counterparty is much less likely to default 
than those names in the portfolio. However, defaults of names in the portfolio increase the 
default intensity of the counterparty. It can be observed that the equity tranche is least 
affected by the counterparty risk while the senior tranche is most affected. This is expected 
as the counterparty is much more likely to default due to the contagion effect from defaults 
of names in equity and mezzanine tranches, and therefore the senior tranche investors require 
higher compensation for increased counterparty risk. 
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4 Conclusions 



In this paper we have suggested a general exponential decay copula contagion mixture model 
which unifies the factor model, copula model, and contagion model. The key advantage is that 
one can investigate the interaction of these models and its pricing impact on basket CDS and 
CDO. The ordered default times can be easily generated with the total hazard construction 
method. We have done some numerical tests and compared basket CDS and CDO rates 
with three copula (product, exponential, and Gaussian) contagion mixture models and found 
that there is no model dominating the others in pricing although one model may affect much 
greatly the pricing of some parts of basket CDS and CDO than the other models do. The 
difference of swap rates computed with different models can be substantial. The exponential 
decay rate has great impact on senior tranche rates and kih default rates for large k. We have 
also compared pricing results when there is contagion counterparty risk. Our conclusion is 
that one has to be cautious in pricing basket CDS and CDO when a particular model is used 
as different models may greatly influence the portfolio loss distribution and can significantly 
affect the resulting swap rates. We should not put all blames on the "misplaced reliance on 
sophisticated maths" for the recent financial crisis, see Cookson (2009). No model is best for 
all purposes. Stress test and other risk control procedures should be in place to withstand 
the potential loss due to the wrong choice of models. 

Acknowledgement. The author thanks Duanpeng Wang for the help in numerical tests 
with C++. The author also thanks the London Mathematical Society for its Scheme 4 grant 
support (Grant 4818). 
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